NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 





DEPARTMENT OF MA&ffiMfttlCAL SCIENCES 
SCHOOL OF SCIENCES MID HEALTH PROFESSIONS 
OLD DOMINION UNIVERSITY 
NORFOLK, VIRGINIA 




f 




PRELIMINARY INVESTIGATIONS INTO THE ACTIVE 

CONTROL OF LARGE SPACE STRUCTURES SOLUTION 

OF THE TIMOSHENKO BEAM EQUATIONS BY THE 
METHOD OF CHARACTERISTICS 

(NASA-CR-163tt08) PRELIMI NARY INVESTIGATIONS N80-29418 

INTO THE ACTIVE CONTROL OF LARGE SPACE 

STRUCTURES: SOLUTION OF THE TIMOSHENKO BEAM 

EQUATIONS BY THE METHOD OF CHARACTERISTICS Unclas 

Final Report, 15 Feb. (Old Dominion Univ., 33/18 28357 

* - ; 
John Tweed 





Q 

d 


Final Report 

For the period February 15, 1976 - August 15, 1980 


Prepared for the 

National Aeronautics and Space Administration I 
Langley Research Center \ 
Hampton, Virginia 


Under 

Research Grant NSG 1279 

Dr. Michael F. Card, Technical Monitor 

Structures and Dynamics Division 







August 1980 



DEPARTMENT OF MATHEMATICAL SCIENCES 
SCHOOL OF SCIENCES AND HEALTH PROFESSIONS 
OLD DOMINION UNIVERSITY 
NORFOLK, VIRGINIA 


PRELIMINARY INVESTIGATIONS INTO THE ACTIVE 

CONTROL OF LARGE SPACE STRUCTURES SOLUTION 

OF THE TIMOSHENKO BEAM EQUATIONS BY THE 
METHOD OF CHARACTERISTICS 


% 

John Tweed 


Final Report 

For the period February 15, 1976 - August 15, 1980 


Prepared for the 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 


Under 

Research Grant NSG 1279 

Dr. Michael F. Card, Technical Monitor 

Structures and Dynamics Division 


eye 


Submitted by the 

Old Dominion University Research Foundation 
P. 0. Box 6369 
Norfolk, Virginia 23508 


August 1980 


TABLE OF CONTENTS 


Page 


ABSTRACT 1 

INTRODUCTION 1 

GOVERNING EQUATIONS 3 

THE ALGORITHM 10 

RESULTS 14 

REFERENCES 16 

APPENDIXES 17 


Figure 


LIST OF FIGURES 


1 


Finite-difference characteristic mesh 


11 


PRELIMINARY INVESTIGATIONS INTO THE ACTIVE CONTROL OF 
LARGE SPACE STRUCTURES -- SOLUTION OF THE TIMOSHENKO 
BEAM EQUATIONS BY THE METHOD OF CHARACTERISTICS 

By 

John Tweed* 

ABSTRACT 

In this report, the method of characteristics is used to develop 
an algorithm for solving the damped, Timoshenko beam equations with 
free-free boundary conditions and prescribed initial data. A number 
of special cases are considered and, where appropriate, comparisons 
are made with known results. 


INTRODUCTION 

The possibility of using large orbiting space structures to help 
meet future needs on Earth has prompted industries and government 
agencies to begin preliminary design studies on a number of such 
structures: for example, large communications antennae are proposed 

as a means of handling the ever-increasing quantities of information 
which are being transmitted around the world. High-resolution soil 
moisture radiometry techniques are being put forward as a method of 
accessing the world's food inventory, and, in an effort to relieve 
the energy shortage, large, orbiting, solar power stations are being 
considered which are designed to transmit the Sun's energy to Earth 
via microwaves. The cost of such structures is enormous and propor- 
tional to the amount of construction material used. Keeping material 
to a minimum will result in a considerable lack of stiffness and per- 
haps in dynamic behavior unlike any structure presently in existence. 
To function properly these structures must be kept free of unwanted 
motions such as those due to the docking of other spacecraft, solar 
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winds, and meteroid showers. An active structural control system 
integrated into the structural design is one method of solving this 
problem. Since the design of such a control system is a difficult 
and sophisticated problem, it is felt that simpler problems should 
be studied with a view to increasing the understanding of the dynamic 
behavior of large structures and giving direction to studies in the 
control of these structures. 

A simple structure which retains many of the important features 
of a large, flexible structure is the Timoshenko beam, which admits 
shear deformation and rotary inertia effects. A previous study (ref. 1) 
shows that, in a Timoshenko beam of infinite length, it is possible 
for the displacement response to an initial disturbance to be a dis- 
persive pulse moving along the beam. This type of response occurs 
if the beam is relatively soft in shear and if the wavelength of the 
disturbance is small compared to the length of the beam. Based on 
studies to date, this type of motion is possible within the range of 
stiffness and mass parameters associated with current proposed con- 
figurations of large space structures. Since the study on beams of 
infinite length was not entirely realistic, two subsequent studies 
(refs. 2, 3) were perfoi'med for beams of finite length. In the first 
of these studies (ref. 2) , the modal method was used to examine the 
dynamic response of an undamped, uniform Timoshenko beam to an initial 
(1 + cos) lateral displacement, which was centrally located on the 
beam. Three types of motion were discerned: bearalike motion, a 

dispersive pulse motion, and wavelike motion. Since only the first 
few modes are necessary to fully understand beamlike motions, such 
motions in large space structures present no new problems and are 
easily controlled in the usual way. However, if the structure does 
not behave like a beam but exhibits the other two types of motion, 
the solution to the control problem is not clear: for example, to 

accurately describe a pulse, many modes are needed. If a controller 
is included for every mode, then an adaptive, active control system 
might be necessary since the frequencies and shapes of that many modes 
(200 or more) cannot be known accurately a priori. Because of the 
difficulty of determining the high modes accurately, it was felt 
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that alternative methods of solving the problem should be examined. 

, Therefore, in the next study (ref. 3) a finite-difference program 

was developed. This program was designed to solve the Timoshenko 
beam equations with distributed loads, velocity damping, free-free 
I boundary conditions and the same initial data as that discussed in 

| reference 2. For large disturbances, the finite-difference solution 

! agreed well with the modal solution but required more CPU time. 

When small wavelike disturbances were considered, a fine resolution 
in the spatial domain was required, and, since the wave speeds were 
fairly large, very small time steps had to be taken. The result 

was a significant increase in CPU time and also in buildup error. 

f' 

Consequently, the finite-difference study was abandoned in favor 
f of the present study, which takes advantage of the hyperbolic nature 

of the Timoshenko beam equations by using the method of 
characteristics. 

GOVERNING EQUATIONS 


The equations governing the behavior of a Timoshenko beam may be 
written in the following form: 
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where 

A * A(x) 
bi » bj(x,t) 
b 2 = b 2 (x,t) 
E * E(x) 

G = G(x) 
g 

I * I(x) 

k g = k g (x) 

L 

M * M(x,t) 
q = q(x,t) 
r = r(x,t) 
s = s(x) 
t 

V * V(x,t) 
w = w(x,t) 

Wi = Wi(X,t) 

w 2 * W 2 (x,t) 

X 

y • Y(x,t) 

G » 0(x,t) 
t ■ <Kx,t) 
o * o(x,t) 


effective cross-sectional area 
transverse damping force/unit mass/unit velocity 
longitudinal damping force/unit mass/unit velocity 
Young's modulus 
shear modulus 

acceleration due to gravity 
second moment of cross-sectional area 
radius of gyration of cross section 
length of beam 
bending moment 

distributed force/unit length 
distributed torque/unit length 
shear correction factor 
time 

shear force 

total deflection 

bending deflection 

shear deflection 

distance measured along span 

shear strain 

total slope 

bending rotation 

mass of beam/unit volume 
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In the case of a free-free beam, we are required to solve equations 
(1) to (6) subject to boundary conditions of the type: 

M(0,t) * V(0,t) - 0 , t > 0 

M(L,t) = V(L,t) * 0 , t > 0 

and to initial conditions such as 

w(x,0) * w°(x) 0 < x < L 

i(/(x,0) = i|»°(x) 0 < x < L 

ff (x,0) * ff(x,0) - 0 , 0 < x < L 

where w°(x) and ip° (x) are prescribed functions of x. 

For computational convenience, we describe the problem in terms of dimension- 
less variables as shown below. 

After choosing a suitable unit of time, t Q , we may introduce a dimension- 
less time variable, t * t/t Q , an< * a dimensionless space variable, 5 » x/L. 
Furthermore, if we let A Q * L 2 , I = L 4 , and assume that A(x) * A Q A(£)» 

I(x) = I 0 I( p )» G(x) = G q G(S), E(x) * E q E( 5) , and p(x) = p q p( 5), we may 
introduce the following dimensionless quantities: 


CU,T) 
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w(CL,Tt Q )/L 

CjU.t) 

s 

w i (CL,rt o )/L , i * 1,2 

♦(«,T) 
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*(5L,Tt 0 ) 
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q(5,x) 
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i ■ 1,2 


^U.t) ■ b i (CL,Tt 0 )t Q , 

C°(0 - w°(CL)/L 

♦°tt) - *°(CL) 



In terms of these dimensionless quantities, the equations (1) through (8) 
take the form 
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with boundary conditions 

M(0,t) * V(0,t) * 0 , t > 0 
M(1,t) * V(1,t) « 0 x > 0 


and 


( 15 ) 
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initial conditions: 


;( C , 0 ) - ;°(0 0 < C < 1 1 

■KC.O) - *°(0 0 < C < 1 

f£<5,0) « f£u,0) = 0 , 0 < 5 < 1 

If we now let i « , and set 

3 T 0 T 



we find that equations (9) through (14) can be written as a linear 
first order system 

* D(0 ~ + B(5,t)u + 1(5, t) , 0 < 5 < 1, t > 0 

with boundary conditions 

M(0,t) « V(0,t) - 0 , t > 0 I 

M(1,t) « V(1,t) - 0 , t > 0 I 

and initial conditions 

£U,0) « *(C,0) * 0 0 < C < 1 

V(5»0) -SAG (C) -*°U)J 0 < 5 < 1 

and 

MU.O) = I I H? (?) 0 < 5 < 1 


( 16 ) 
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where 



( 21 ) 



( 22 ) 
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Since the D matrix has eigenvalues 
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and corresponding eigenvectors 
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(cont'd.) 
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(25) 

(concl'd. ) 


which are linearly independent, the system is seen to be hyperbolic. On 
solving this hyperbolic system for the vector u we calculate c(?,t) and 
$(C,t) from th* eouations 


C(C.t) 


» ;°(C) +£ C(c.y) dy 


( 26 ) 


and 


♦ (5 .t) - *°(C) *jT T Hi.y) dy 


(27) 


The matrix T whose columns are the eigenvectors xj, x 2 , x 3 , and x 4 

of D is nonsingular and such that 

T* 1 DT * A * diagCXj.Xj,^,^) (28) 

Therefore, if we introduce the new variable v * T 1 u and let 

B* • T _1 BT and f* * T _1 ?, we find that equation (10) may be written in 

the form 


3 v 
3t 


A(0 ♦ B*(C,t)v ♦ ?*(C.t) 
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while the boundary conditions become 


Vi(0,t) ♦ v 2 (0,t) ■ 0 
v 3 (0,t) + v 4 (0,r) - 0 
Vi(1,t) ♦ v 2 (1,t) « 0 
v 3 (l,t ) ♦ v 4 C1,t) « 0 


t > 0 


and the initial conditions become 


v(£,0) > v°(£) 


SAG 

r 0 
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0 < £ < 1 


THE ALGORITHM 


On the characteristic curves C. , which are defined by the 
differential equations ^ - -X., i * 1,2, 3, 4, equation (30) takes the 
simple form 


dv . 

l 



+ 


f i ) j7 


Equation (33) forms the basis of an algorithm for solving the system 
(30) numerically. This algorithm, which has been discussed by Leonard 
and Budianskv (refs. 4,5) for systems without damping, is described 
below in the case of a uniform, homogeneous beam with square section. 
In this case the eigenvalues reduce to ±aj/ /s , +a 2 . and hence thi 
characteristics become straight lines. 




On the strip 0 < 5 < 1, t > 0, we construct a mesh of size A£, 

At where At ■ AC/a2* We let P * <*\/<* 2 A and make use th® notation 
5^ « iA£, t ^ a jAx, f(C i ,Tp ■ f?, eve. A typical cell in this mesh 
is shown in figure 1. If, for a fixed value of j, the solution 
is known at the mesh points (5^»Tj) f it can be found at the points 
Uj + p,^) by using three-point Lagrange interpolation. The solution 
may then be advanced in time to the points (5^ + 1) by means of 

equation (341 below. 



Figure 1. Finite-difference characteristic mesh. 
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which is a finite-difference expression for equation (33) with 

taking the values p, -p, 1, -1 as ra takes the values 1, 2, 3, 
or 4, respectively. Subject to the boundary conditions 

( V 1 + V 2 )j +1 * ( V 3 + V 4)o +1 " 0 


+ vJV + 1 = fv, + v.V + 1 * 0 


j > 0 


( v l + V 2)i +1 = ( V 3 T * k J n 
the solution of equation (34) may be written in the form 

?|*> * cl* 1 ?|* 1 . v° - o < i < », j > o 

_ -H+i . . . 

where the mth component of r< is given by 
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if 0 < i < n, by 
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if i = n, where in all cases 

2 — 

Atr j + l u Axai AtsAcx 2 

3 ■ T b >i -Jr • 

d = ^b 2 | + ^, h = 1 + 2a, and 
k = (1 ♦ 2d) + 4bc{ 0 i 6 ni' 


The quantities of physical interest are obtained from the solution of 
equation (36) with the help of the formulas 







-v 2 )^ +1 ♦ (V 1 - v,)i] 

- v 4 )^ +1 + (v 3 - voj] 


(39) 


(40) 


(41) 
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(43) 



(Vl ♦ v 2 )| +1 


and 



(v 3 
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RESULTS 


The algorithm described above has been incorporated into a 
computer code for solving the damped, Timoshenko beam equations with 
free- free boundary conditions and prescribed initial data. A number 
of special cases have been investigated and, where appropriate, 
comparisons have been made with earlier results. In order to facilitate 
these comparisons, the time scale t was taken to be the period 
of the first mode of an equivalent Euler beam (ref. 7) : 

\h 


l o - 


2rrL" 


(4 .730) 2 [ E o I 


(45) 


and results were obtained for different slenderness ratios, SLR = 

h h 

L(A/I) , and velocity ratios, V D = (G /E ) . 

K 0 0 

The graphs in appendixes A through F show the undamped displace- 


ment response tj initial data of the form 

c(5,0) = jjl + cos {x U " 3>}] H [* 2 - ^ - U 2 ] (46) 

♦ (5,0) = - y sinjyU - J)}h£* - (2g - 1) 2 J (47) 

V(5,0) = 0 (48) 

and 

M(5,0) = 0 (49) 

where H(x) =jl, x>o) is the heavyside function. 

(0, x<0) 
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Appendixes A and B deal with the cases V R ■ 0.62, SLR * 50.0, and 
A * 0.5 and 0.1, respectively. The results given here agree with those 
obtained by the modal method (ref. 2) and also by the direct i'inite- 
difference method (ref. 3) . A comparison of machine time used shows 
that the method of characteristics is about 20 times faster than the 
direct finite-difference method. It appears, however, that the 
characteristic method is not as efficient as the modal method 
especially when short wave lengths are involved. This is because 
the high spatial resolution required to see short waves demands that 
very small time steps be used, and the resulting buildup error quickly 
swamps out the desired solution. Appendixes C, D, E, and F show 
corresponding results for the cases V R = 0.3, SLR = 50.0, X = 0.5; 

V R = 0.3, SLR = 50.0, A = 0.1; V R = 0.1, SLR = 50.0, A = 0.5, and 

V R = 0.1, SLR = 50.0, A = 0.1, respectively. 

Appendixes G through K show how the displacement, bending angle, 
shear and bending moment respond to initial data of the type represented 
by equations (46) to (49) in the presence of uniform velocity damping. 
Results are presented for the cases bj = 0.0, b 2 = 0.0; bj = 1.5, 
b 2 = 0.0; bi = 3.0, b 2 = 0.0; bj = 0.0, b 2 = 3.0, and bj = 0.0, 
b 2 = 12.0, respectively, with V D = 0.62, SLR = 50.0, and A = 0.2. 

A comparison of the solutions obtained shows that transverse damping 
is much more effective than rotational damping. 

An attempt was made to analyze local damping and point- load-rate 

effects. While results due to Leonard (ref. 5) could be reproduced 

approximately, in general the method proved to be poor for this class 
of problems. Again it is found that the high spatial resolution 
requires results in very small time steps and buildup error which 
quickly swamps out the desired solution. 
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APPENDIX A 


Displacement response to a (1 + cos) initial disturbance; 

V_ * 0.62, SLR - 50, X = 0.5. 
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APPENDIX B 

Displacement response to a (1 ♦ cos) initial disturbance; 
V R = 0.62, SLR = 50, X = 0.1. 
















APPENDIX C 


Displacement response to a (1 + cos) initial disturbance; 


V_ = 0.3, SLR = 50, \ = 0.5 
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APPENDIX E 

Displacement response to a (1 + cos) initial disturbance; 

V D = 0.1, SLR = 50, X = 0.5. 
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APPENDIX F 


Displacement response to a (1 ♦ cos) initial disturbance; 
V_ * 0.1, SLR = SO, \ - 0.1. 
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APPENDIX G 


Response to a (1 + cos) initial displacement; 

V D = 0.62, SLR = 50, A = 0.2, bj = 0.0, b 2 = 0.0. 
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APPENDIX H 

Response to a (1 + cos) initial displacement; 

V D = 0.62, SLR = 50, A = 0.2, b x = 1.5, b 2 = 0. 
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APPENDIX I 

Response to a (1 ♦ cos) initial displacement; 

V R = 0.62. SLR = 50, \ = 0.2, bj = 3.0, b 2 = 0.0. 
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APPENDIX J 


Response to a (1 ♦ cos) ’ itial displacement; 

* 0.62, SLR = 50, X = 0.2, bj * 0.0, b2 = 3.0. 
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APPENDIX K 


Response to a (1 + cos) initial displacement; 

V_ = 0.62, SLR = 50, X = 0.2, b! = 0.0, b 2 = 12.0. 
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